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ABSTRACT 

We present results from a low resolution (R ~ 300) near-infrared spectroscopic variability survey 
of actively accreting T Tauri stars (TTS) in the Taurus-Auriga star forming region. Paschen and 
Brackett series H i recombination lines were detected in 73 spectra of 15 classical T Tauri systems. 
The values of the Paniip/Pa/3, Bmiip/Br7, and Br7/Paniip H i line ratios for all observations exhibit a 
scatter of < 20% about the weighted mean, not only from source to source, but also for epoch-to-epoch 
variations in the same source. A representative or 'global' value was determined for each ratio in both 
the Paschen and Brackett series as well as the Br7/Pant,p line ratios. A comparison of observed line 
ratio values was made to those predicted by the temperature and electron density dependent models 
of Case B hydrogen recombination line theory. The measured line ratios are statistically well-fit by a 
tightly constrained range of temperatures (T < 2000 K) and electron densities (10^ < no < 10^° cm~'^). 
A comparison of the observed line ratio values to the values predicted by the optically thick and thin 
local thermodynamic equilibrium cases rules out these conditions for the emitting H i gas. Therefore, 
the emission is consistent with having an origin in a non-LTE recombining gas. While the range of 
electron densities is consistent with the gas densities predicted by existing magnetospheric accretion 
models, the temperature range constrained by the Case B comparison is considerably lower than that 
expected for accreting gas. The cooler gas temperatures will require a non-thermal excitation process 
(e.g., coronal/accretion-related X-rays and UV photons) to power the observed line emission. 

Subject headings: circumstellar matter — infrared: stars 



1. INTRODUCTION 

In the decades since their discovery, photometric vari- 
ability continues to be a defining characteristic of T 
Tauri stars (TTS) and has been observed at most wave- 
lengths with timesc ales as short as hours (e.g.. lJovlll945l : 
iBouvier et al.lfl993l ). In addition to the widely studied 
photometric variations of TTS, astronomers have begun, 
more recently, to probe the variability of the unique spec- 
tral emission features associated with these sources us- 
ing multi-epoch spectroscopic observations to note vari- 
ations in the s hapes and line strengths of the emission 
features (e.g., iJohns fc Basril Il995al: iBasri et al.l Il997t 
lAlencar et al.l 120011 : IBouvier et al.ll2007t ). An investiga- 
tion of the physical conditions that give rise to these 
emission features will illuminate many of the processes 
involved in the formation of these stars. 

Single-epoch observations of H i Balmer series line 
emission from TTS reveal a complicated kinematic line 
structure, with features often including blue-shifted 
and/or red-shifted absorption components. As evi- 
dent from the catalog of Ha features presented in 
iReipurth et al.l (119961 ). a variety of complex profiles ex- 
ist for this feature and several other Balmer emission 
lines. The same spectral line studies reveal quite large 
full-width at half of maxima (FWHM) for the emission 
components of these features and suggest that the gas 
is moving at hundreds of km s~^ relative to the sources. 
The complexity of these emission features alone has led 



Electronic address: jsb3r@virginia.edu, spm5x@virginia.edu, mfs4n(t 
^ Department of Astronomy, University of Virginia, P.O. Box 

400325 Charlottesville, VA 22904-4325 

^ Harvard-Smithsonian, Center for Astrophysics, 60 Garden St., 

Cambridge, MA 02138 



stars: pre-mam sequence 



to different conclusions about the origins of this emit- 
ting gas. However, the blue-shifted absorption compo- 
nents or P Cygni profiles for these emission features are 
most likely produced by gas entrained in an outfl owing 
wind ("Hartmann ct al.l(l990l : iGrinin fc Mit skevichI fl99ll : 
[Calvet ct al. 1992). As a result, the observation of a large 
Ha equivalent width (EW) for an individual source was 
initially interpreted as evidence of and arising entirely 
from an active outflowing wind. 

In another model of the star-disk interaction, the H i 
lines and other TTS phenomena were formed in a viscous 
boundary layer, w here the accretion disk touched the 
surface of the star (jLvnden-Bell fc Pringlel[l97l : iBertoutJ 
Il989| ). However, the inability of both the boundary layer 
and stellar wind models to explain the shapes and widths 
of H I emission features led to the ap plication of the 
magnetosph eric accretion model to TTS (jGhosh fc LambI 
119791 : lKonigI.1991.) . Previously developed to explain disk 
accretion onto compact objects, the magnetospheric ac- 
cretion model established a new paradigm for the star- 
disk interaction in these pre-main sequence stars. This 
model involves a stellar magnetosphere that truncates 
the disk, opening a gap with a size of a few stellar radii 
between the disk and the surface of the star. Loop- 
ing stellar magnetic fleld lines thread the inner most re- 
gions of these disks, acting as accretion funnels through 
which disk matter may free fall onto the stellar sur- 
face. The free falling material, confined to a magnetic 
channel, naturally explains the large velocities observed 

ijfara.Ehfl, jBwSz«5i§aMte(fe,aapea-^oae(»geter^arda8u,agpea avirginia.edu] 

ance of red-shifted absorption components. In addition, 
the accretion shocks formed at the point of impact at 
the stellar surface provide the heated gas necessary to 
produce the 'blue' continuum excesses characteristic of 
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actively accreting TTS. Several authors have employed 
the magnetospheric accretion model to predict the fluxes 
and shapes of a few of the op tical an d infrared H i 
emissi on features (jHartmann et a l. 1994; MuzeroUe et al.1 
Il998airbl ). These modeling efforts support the idea that 
the bulk of the emitting H i gas is confined to the columns 
of accreting gas, rather than an outflowing wind. 

Advances in infrared spectroscopy have allowed as- 
tronomers to observe at high spectral resolution a num- 
ber of H I Paschen and Balmer series emission lines, 
predominantly Pa/3 an d Br'- y, from TTS (Naiita et al.l 
[1996; Muzcrollc et al. 1998al |bl: iFolha fc Emersonlll998l . 
[2001. : iNatta et al.ii2006i) . These infrared H i hues, which 
possess lower opacities relative to the Balmer series, 
'conspicuously' lack the blue-shifted absorption features 
present in Ha and H/3 profiles (Folha fc Emerson 1998) . 
The infrared lines also appear to correlate with the 
mass accretion rates determ ined from measurements of 
the 'blue' continuum excess (|Gullbring et al.lll998L[2000l : 
iMuzeroUe et al.ll2001f ) suggesting that the H i emission 
is related to accretion activity. Therefore, as a result of 
the successes of the magnetospheric accretion model in 
explaining the shapes of the H i emission features and 
the correlation of these features with other accretion di- 
agnostics, the interpretation of the H i emission lines and 
the location of the emitting gas has shifted in the last 15 
years. Where once the presence of a strong Ha emission 
line in the spectrum of a TTS was widely interpreted as 
an indicator of an outflowing wind and the 'blue' con- 
tinuum excesses were evidence of a hot boundary layer, 
today the H i features and the 'blue' excesses are both 
interpreted as evidence of magnetospheric accretion. 

More recently, spectro-astrometric observations of sev- 
eral classical TTS (cTTS; TTS possessing evidence for 
actively accreting circumstellar disks) and Herbig AeBe 
stars possessing symmetric (Type I) Pa/3 emission fea- 
tures have shown that some of the emission in the line- 
wings have measurable offsets with respect to the cen- 
tral sources. These results suggest that outflowing gas 
is a li kely sou rce fo r this component of the emission 
(jWhela n et al. . l2004l) a r id po tentially explain the diffi- 
culty |Fbllia^.^^merso3 (|2001h encountered modeling the 
line-wings of these emission features with a magneto- 
spheric accretion scenario that did not include a wind. In 
addition, AMBER/VLTI observations of Herbig Ae star 
HD 104237 revealed that the Br7 and _K-ba nd continuum 
emiss io n have simila r spati al distributions ("Tatu Ui et al.l 
|2007|) . iTatulli et all (|2007f l successfully explain the vis- 
ibility of the Br7 emission by modeling it with an out- 
flowing wind and places the emission at the base of the 
jet with a radial extent of 0.2 to 0.5 AU. Observations 
such as these indicate that a portion of the near-infrared 
H I emission features can be produced in outflowing gas. 

Despite these apparent advancements in understand- 
ing the origin of the H i emission features, the physi- 
cal conditions of the emitting gas (e.g., ionization state, 
level popula tions, t emperatur e, de nsity) remain poorly 
constrained (lAlencar fc Basril [20001: Iwhelan et al. 2004 : 



iTatufli et al] 120071 : lEisneill2007D . iMartinI (|l99d rcalcu- 
lated a temperature and density profile for gas traveling 
along an accretion funnel from the disk to the star by 
modeling several sources of gas heating (e.g., ambipolar 
diffusion, adiabatic compression, photoionization, and 
external radiative excitation) and cooling (e.g., Ca ii 



and Mg II ions). Assuming an initial gas temperature 
of 3000 K for gas located at the inner edge of the disk 
where it is loaded into the accretion funnel, for a fldu- 
cial m ass accretion rate of M = 10~^ Mq yr~^, iMartinI 
(1996) found that the gas temperature proflle peaks at 
T - 6500 K. 

In another investigation seeking to constrain the 

temperature of magnetospherically accreting gas, 

[MuzeroUe et all ([2001) calculated Ha and H/3 line 
fluxes and profiles dependent upon Al and a range 
of gas temperatures. The authors found 'optimal' 
temperatures of 6000 < T < 12000 K for a range of 
accretion rates 10"'' > M > 10"^" Mq yr-\ with the 
lower gas temperatures corresponding to sources with 
higher mass accretion rates. This s tudy su g gested that 
the gas temperatures predicted by IMartinI (|1996[ ) were 
too low to produce the observed line luminosities for the 
M values considered. These temperatures agree with 
the standard thermal gas temperature of T ~ 10^ K 
observed for most H II emission regions. 

A direct measurement of the physical conditions of the 
gas will provide a valuable test of the models and may 
help to further constrain the origin of the H i line emis- 
sion. In optically thin gases, line ratios provide valu- 
able information about the physical conditions. How- 
ever, at high densities, optical depth effects require ra- 
diative transfer modeling and thus complicate the inter- 
pretation of these line measurements. In cTTS, Ha emis- 
sion lines often appear to be optically thick and may be 
formed by the contributions of emitting gas in multi- 
ple regions of the star-disk syst em, not lim ited to the 
accretion flow alone (Kurosawa et al.l |2006[ ) . The near- 
infrared H I lines, with higher n-values, have lower op- 
tical depths and ap pear to be less cornplicat ed than the 
Balmer transitions (jFolha fc Emersonll2001h . Observa- 
tional attempts to compare the measured values of the 
Pa/3/Br7 H i infrared line ratio to the values predicted 
for optically thin and thick hydrogen gas in LTE, in 
the past, hav e failed to determine a temperature for the 
emitting gas ([MuzeroUe et al.|[2001h . Similarly, compar- 
isons of the Pa/3/Br7 ratio to the values predicted us- 
ing the Case B approximation of hydrogen recombina- 
tion line theory ([Storev fc Humnieii (|1995af ) and refer- 
ences therein) have also failed to flnd a meaningful re- 
sult. In fact, the lack of agreement between the observed 
and predicted line ratios lead to the conclusion that the 
emitting H i gas is neither optically t hin or thick nor 
can it be described as a recombining gas 'MuzeroUe et all 
(|2001[ ). However, in most cases where comparisons were 
made, the values for this line ratio were determined 
from spectra collected at different times. Given that 
the variability of i nfrared H i line fluxes has timescales 
as short as hours (iJohns fc Basril [l995bl : lAlencar et al.l 
1200 It iBouvier et al.ll2007( ). a Une ratio calculated from 
two distinct, non-simultaneous spectral observations will 
likely be misleading. Therefore, it is not surprising that 
previous comparisons of a single line ratio have failed 
to meaningfully constrain the physical conditions of the 
emitting gas. S imultaneous obser vations of not only Pa/3 
and Br7 (e.s... iGatti et al.ll2006D . but including several 
other Paschen and Brackett series emission lines will im- 
prove upon these previous comparisons. 

We present an analysis of H i spectral line data simul- 
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TABLE 1 
Target Star Data 



Star Name 



Ha EW 

(A) 



M 



Lx' 



(W'^Mq yr-i) (lO^O erg s"!) 



AA Tau 


37-46 


0.04 


1.039 


BP Tau 


40-92 


0.23 


1.482 


CWTau 


135-140 


1.27 




CYTau 


55-70 


0.08 


0.108, 0.202 


DP Tau 




2.08 




DG Tau 


63-125 


9.19 


0.252 


DM Tau 


139 




0.181 


DO Tau 


100 






DQ Tau 








DRTau 


106 


19.01 




PP Tau 


38 






GK Tau 


15-35 


0.14 


1.244 


HLTau 


43-55 




3.285 


RYTau 


13 




5.242 


TTau 


41-60 




9.395 


UY Aur 









'^ Compiled fron i literature by White & Glicz (2001). *" Val- 
ues measured inlCalvet fc GuUbrina (.1998. ) and compiled by 

[Johiis-Krull & Gafford (2 0021 ).'^ Compiled from literature by 

IGiidel fc Telleschi (2007) . 

taneously acquired as part of a multi-epoch near-infrared 
spectroscopic variability study of cTTS in the Taurus- 
Auriga star forming region. Multiple Paschen and Brack- 
ett emission features were detected in the spectra of our 
sample of actively accreting cTTS collected during the 
three year survey. We determine the values for 16 distinct 
line ratios that are representative of the whole sample of 
TTS. In SjHl we compare the set of observed ratios, broken 
into three groups, Paschen series, Brackett series, and 
Br7/Pan„p, where Uup is the upper level of the transi- 
tion, with the temperature and electron density sensitive 
Case B models as well as the temperature-only sensitive 
optically thick and thin LTE scenarios. The results of 
these comparisons and their predictions of the physical 
conditions of the emitting H I gas are discussed in SJH 

2. OBSERVATIONS 

Multi-epoch spectroscopic observations were made 
of 16 cTTS in the Taurus-Auriga star forming re- 
gion. The spectra were acquired using CorMass 
(Cornell-Massachusetts), a cross-dispersed low resolu- 
tion (R ~ 300) spe ctrometer wit h continuous coverage 
from 0.8 to 2.5 ^m (Wilson ct al. 2001), with the 1.8-m 
Vatican Advanced Technology Telescope (VATT) atop 
Mt. Graham in Arizona. The wavelength coverage of 
CorMass permits simultaneous observation of multiple 
Paschen and Brackett series H i emission lines as well as 
He I (1.083 ^m) and the Ca ii infrared triplet (0.8498, 
0.8542, and 0.8662 ^m). A total of 104 spectra of 16 tar- 
gets were obtained during the period of December 2003 
to January 2005. The cTTSs observed, along with pre- 
viously measured values for Ha EWs, M, and X-ray lu- 
minosities are listed in Table 1. 

On VATT, CorMass slit dimensions are 1'.'6 in width 
and 12'.'2 in length. Observations of both the target 
cTTSs and telluric calibration sources were made us- 
ing a standard ABBA nod pattern, in which the source 
was centered in the slit at position A, observed, then 
'nodded' or moved along the slit to position B, observed 
twice with exposure times equal to the first exposure. 



then nodded back to A for the fourth and final expo- 
sure. The nod pattern allows for efficient removal of 
thermal background, sky emission lines, and dark cur- 
rent. Data reduction was performed using the IDL pro- 
gram CORMASSTOOL adapted froin the Spextool re- 
duction package ( Gushing et al.l[2004[ ). Bad pixel cor- 
rection, flat fielding, and wavelength calibration are ac- 
comphshed through the use of the XCORMASSTOOL 
routine. The wavelength calibration package utilizes re- 
solved OH telluric emission features to individually cal- 
ibrate each spectrum. G2V main sequence stars (HD 
26710, HD 283697, HD 56513, and HD 30455) were used 
as telluric calibration sources. All H i Paschen and 
Brackett series absorption lines not blended with tel- 
luric features were removed from the calibrator spectra 
by linear interpolation prior to performing the telluric 
correction. The error introduced from the interpolation 
procedure is estimated to be less than 1% of the contin- 
uum level. The corrected target spectra were multiplied 
by a blackbody function matching the temperature of 
the calibration sources in order to preserve the shape of 
the spectral energy distribution. Figure 1 contains seven 
multi-epoch observations of DR Tau, spanning the dura- 
tion of the survey carried out with CorMass at VATT. 

The H I line fiuxes were determined using a standard 
IRAF Gaussian fits procedure. T he fluxes w ere then 
dereddened using Ay values froi n I Strom et a l. (1989. 
iKenvon fc Hartmannl (|1995[) . and IWhite fcGhea (,200 ll 
and a standard A>, ^ A~^ -^'* extinction law adopted from 
iMartin fc Whittetl ([l99l . 




Pig. 1. — As an example of the data collected for the multi-epoch 
survey, seven of the spectra obtained for the source DR Tau are 
presented above. Each spectrum contains multiple H I recombina- 
tion emission features, the Ca II infrared triplet in emission, and 
the He I (1.083 fira) feature, seen here in absorption. Individual 
spectra are labeled with the UT date on which it was observed. 
Vertical dotted lines locate the spectral features. 



3. ANALYSIS OP H I LINE RATIOS 

The line fluxes for nine Paschen (Pa/3, 7, 6, 8, 9, 10, 
11, 12, and 14) and eight Brackett series (Br7, 10, 11, 12, 
13, 14, 15, and 16) lines were measured. These H I line 
fluxes^ appear to differ by greater than an order of mag- 

^ The large width of the CorMass slit would enable us to very 
accurately flux calibrate the spectra under photometric observing 
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nitude from source to source. Multi-epoch observations 
of individual sources reveal significant line variability as 
can be determined by inspection in Figure 1. 

While the magnitudes of the line fluxes and their vari- 
ability are sensitive to uncertainties in the flux calibra- 
tion, the line ratios are independent of such errors. In 
regimes where the gas is either optically thick or thin in 
LTE or is well-described by the approximations of hydro- 
gen recombination line theory, these ratios can provide 
an accurate determination of the physical conditions of 
the emitting gas. The spectra in this survey have the 
advantage over previous studies of H i line ratios ob- 
served for TTS because they allow for a large number 
of emission lines to be measured simultaneously. Here 
we describe the determination and present an analysis 
of three different series of line ratios and compare these 
observed values to those predicted by three distinct ra- 
diative regimes in an attempt to constrain the physical 
conditions of the emitting gas. 

3.1. Calculating H I Line Ratios 

Figure 2 shows the dereddened line fluxes, including 
observational uncertainties, for Pan„p to Pa/3 and Brri„p 
to Br7. Only emission lines measured to be 3a detec- 
tions or better are plotted and only these line are used 
in this analysis. Seventy-three of the 104 spectra have at 
least one useful ratio, and only lines from these spectra 
are used in this analysis. None of the observations of RY 
Tau have 3a detections of two H l lines, but all other stars 
in our sample have at least one spectrum in the set of 73. 
As shown in Figures 2, a strong linear correlation exists 
between the H i line fluxes with surprisingly little scat- 
ter, considering that the data represent many different 
sources and multiple epochs of observations of the same 
source. In Figure 3, Pa/3 equivalent widths are plotted 
versus the values of the Pa7/Pa/3 line flux ratios illustrat- 
ing the intra-source variability and the consistency of the 
values measured for the line ratios. Table 2 provides the 
values for the equivalent widths and flux ratios for the 
sources designated by colored points in Figure 3 along 
with the Julian date for each epoch of observation. For 
sources such as CY Tau, BP Tau, and DR Tau that pos- 
sess large epoch-to-epoch variations in the values of their 
line fluxes and substantially different average flux values, 
the values of the line ratios remain constant within the 
uncertainties. A similarly small scatter in the values of 
this line ratio is observed for sources such as T Tau and 
DG Tau that possess relatively small variations in the 
magnitudes of their line fluxes. 

The relatively small scatter in the observed values of 
the line flux ratios for this sample of 15 T Tauri systems 
possessing a broad range of intrinsic mass accretion rates 
(see Table 1) suggests that the ratios are not sensitive to 
variations in M. For this reason, and as a first step to- 
ward measuring the average physical conditions of the 
gas in these systems, we determine a 'global' line ratio 
representative for all of the stars in this survey for each 
ratio included in the analysis. In order to determine the 
representative line ratios, an average value for each ratio 

conditions. However, some observations were made during inferior 
conditions introducing error into our flux calibration and line flux 
measurements. This discussion of line flux variability refers to the 
large scale variations attributed to source variability and not the 
uncertainties in the flux calibration. 
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Fig. 2 — For each spectra in our sample for which we detect 
these Paschen series lines, we plot the measured Pa/3 line flux ver- 



ated measurement uncertainties and find a strong linear correla- 
tion for all emission features. Both axes of all plots are in units of 
10^^^ ergs cm~^ s~^ . Each solid line plotted through the data has 
a slope equal to the weighted average ratio of the data points and 
represents the 'global' line ratio for each plot. The dashed lines 
represent the estimated variance and provide a useful measure of 
the scatter in each of the line ratios shown on each graph. Also 
included are similar plots for the Brackett series lines, using Br7 
and seven other Brn„p transitions. 



weighted by the uncertainty was found. In addition, the 
statistical "estimated variance," was calculated to quan- 
tify the scatter about the average line ratio values. In 
Figure 2, the solid lines have slopes given by the 'global' 
line ratios determined for each pair of line fluxes. The 
dashed lines have slopes equal to 'global' line ratio plus 
and minus the estimated variance and thus represent the 
scatter in the measurements. 

Each ratio uses the Paschen or Brackett series transi- 
tion with the lowest excited state as the reference tran- 
sition (Pa/3 and Br7). The number of stars with de- 
tectable H I features diminished for the transitions with 
higher n-values. As a result, the representative line ra- 
tios involving the highest transitions in the Paschen and 
especially the Brackett series are determined by only a 
few of the TTS in our sample. However, the small intrin- 
sic source-to-source scatter observed for the other ratios 
that included larger samples of line fluxes suggest that, 
even in the cases where only a few sources are used to 
determine the value of the line ratio, it is reasonable to 
infer that the ratios are representative of the sample as a 
whole, within the uncertainties. We proceed by making a 
direct comparison between the observed values for these 
sets of line ratios and the predictions of line formation 
theory. 
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TABLE 2 
Selected line data. 
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Fig. 2 — Continued 
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Fig. 2. — For each of the 73 spectra with detectable Pa/3 and 
Pa7 emission features, we plot the values of the Pa/3 equivalent 
widths versus the values of the Pa7/Pa/3 line flux ratios. Sources 
for which we have six or more different epochs of observations have 
been color coded. The dark horizontal line in the center represents 
the weighted average of the line ratio for all of the data points. The 
dashed horizontal lines represent the "estimated variance" of the 
line ratio value. For large variations in the magnitude of the Pa/3 
EWs, relatively small variations are seen in the values of the line 
ratio. The colors highlight this result for the individual sources. 



3.2. Comparison to Case B models 
3.2.1. Case B 

The determination of the electron cascade and the pop- 
ulation densities of the energy levels in recombining hy- 
drogenic atoms is a well-studied problem that, histori- 
cal ly, has been divided in to two distinct Cases, A and 
B (jBaker fc Menzellll938[ ). Case A applies to a gas for 



Star 



Epoch 
(JD-2450000) 



Pa/3 
EW{A) 



Pa7/Pa/3 
flux ratio 



BPTau 


2976.8 


-5.67 ± 0.91 


1.00 ± 0.36 




2977.8 


-6.86 ± 0.83 


1.01 ± 0.35 




2986.7 


-8.37 ± 0.73 


1.11 ± 0.21 




2986.9 


-9.2 ± 1.1 


0.77 ± 0.20 




2987.7 


-6.54 ± 0.84 


0.96 ± 0.26 




2988.7 


-9.77 ± 0.78 


0.93 ± 0.16 




2989.6 


-8.1 ± 1.1 


0.96 ± 0.28 




2989.9 


-10.6 ± 1.6 


0.81 ± 0.26 




3065.7 


-13.0 ± 1.1 


0.92 ± 0.21 




3074.7 


-5.84 ± 0.82 


0.81 ± 0.22 




3301.8 


-14.8 ± 1.1 


0.78 ± 0.12 




3399.7 


-13.96 ± 0.97 


0.94 ± 0.13 




3401.8 


-15.64 ± 0.79 


0.897 ± 0.096 


CWTau 


2978.9 


-13.02 ± 0.89 


0.82 ± 0.16 




2986.7 


-12.04 ± 0.63 


0.716 ± 0.089 




2987.6 


-12.58 ± 0.89 


0.81 ± 0.15 




2988.7 


-13.67 ± 0.47 


0.713 ± 0.067 




2989.6 


-16.1 ± 1.3 


0.75 ± 0.14 




3301.8 


-17.39 ± 0.53 


0.654 ± 0.045 




3399.7 


-12.86 ± 0.86 


0.85 ± 0.14 


CYTau 


2978.8 


-4.2 ± 1.3 


0.98 ± 0.46 




2988.7 


-3.77 ± 0.68 


0.99 ± 0.37 




2989.6 


-11.3 ± 1.0 


1.06 ± 0.19 




3074.7 


-5.8 ± 1.4 


0.86 ± 0.35 




3301.8 


-6.60 ± 0.82 


1.09 ± 0.28 




3399.7 


-4.2 ± 1.1 


2.1 ± 1.0 


DGTau 


2978.9 


-12.71 ± 0.86 


1.02 ± 0.16 




2986.9 


-13.69 ± 0.58 


1.011 ± 0.096 




2987.9 


-15.96 ± 0.79 


0.81 ± 0.10 




2988.7 


-12.15 ± 0.54 


0.931 ± 0.099 




2989.7 


-9.69 ± 0.69 


0.96 ± 0.16 




3301.9 


-12.58 ± 0.82 


0.88 ± 0.13 




3399.8 


-12.5 ± 1.4 


0.91 ± 0.19 


DRTau 


2978.9 


-23.37 ± 0.67 


0.909 ± 0.053 




2986.8 


-23.64 ± 0.34 


0.889 ± 0.027 




2986.9 


-23.05 ± 0.77 


0.884 ± 0.059 




2987.7 


-22.86 ± 0.89 


0.896 ± 0.070 




2988.7 


-17.20 ± 0.52 


0.802 ± 0.054 




2989.7 


-15.4 ± 1.8 


0.79 ± 0.19 




2989.9 


-16.5 ± 1.2 


0.83 ± 0.13 




3074.7 


-21.89 ± 0.45 


0.854 ± 0.039 




3399.6 


-19.08 ± 0.64 


0.881 ± 0.061 




3401.8 


-18.7 ± 2.8 


0.95 ± 0.24 


TTau 


2979.9 


-7.44 ± 0.37 


0.92 ± 0.12 




2986.9 


-7.63 ± 0.73 


1.03 ± 0.22 




2987.9 


-11.6 ± 1.1 


0.74 ± 0.14 




2988.7 


-9.39 ± 0.68 


0.89 ± 0.15 




2989.7 


-9.4 ± 1.5 


0.89 ± 0.33 




3301.9 


-10.36 ± 0.75 


1.01 ± 0.17 




3401.8 


-8.81 ± 0.88 


0.82 ± 0.15 



which every transition in the hydrogen atom is optically 
thin such that any photons emitted by atoms within 
the gas escape from the emitting region with no sub- 
sequent scattering even ts. The Case B approximation of 
iBaker fc Menzell ()1938[ ) assumed the gas to be optically 
thick to the ultraviolet photons associated with the Ly- 
man series and continuum and optically thin to photons 
from all other transitions. These early models included 
only radiative capture and cascade and, as a result, were 
only sensitive to the gas temperature. In the decades 
since this work, improvements upon their estimates of 
the n, /-level populations came with the inclusion of 
collisional effects (i.e., collisional ionizat ion, thr ee-body 
recombinatio n , and all A n transitions; ISeatonI ([1964); 
iBrocklehurstI (jl970l Il971h ). The introduction of col- 
lisional effects adds an electron density sensitivity to 
the relative n, Z-level populations of the hydrogen atoms. 
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The refore, the current stat e-of-the-art Case B mod- 
els (jHummer fc StorevI [l987t IStorev fc Hummed Il995a^ 
compute the relative intensities for H I recombination 
lines for a range of temperatures (500 < T < 30000 K) 
and electron densities (100 < no < 10"'^^ cm~^). 

Case A is applicable to only the lowest density nebula 
and astrophysical environments and, most likely, will not 
be the most appropriate approximation for gas located 
in the very inner regions of an accreting T Tauri system 
(| OsterbroclJ 11983) . Therefore, in the following analysis, 
we focus on the Case B approximation. 

3.2.2. Panup/Pa/3, Urnup/BvY, and Hr^f/P&riup 



A Fortran progra m 

IStorev fc Hummeij (|l995allbl ) 

nied by the necessary data 



provided by 

and accompa- 

fi-les (on-line at 



http://vizier.u-strasbg.fr /viz-bin/VizieR?-source=VI/64 1 
was used to produce models ol the Paschen and Braclcett 
decrements as well as the line ratios of Br7 to Pan„p 
for a grid of 130 different temperature and electron 
density combinations for 500 < T < 30,000 K and 
100 < Uo < 10^"* cm^"^. For a statistical comparison 
between the Case B models and each series of observed 
line ratios, the reduced Chi-square (x^) value for each 
model was computed. 
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TABLE 3 
Case B Models with Smallest x^ 



Line Ratio 


nc 


T 


Reduced 




(cm-3) 


(K) 


x' 


Panup/Pa/3 


IQlO 


1000 


0.93 




IQlO 


3000 


0.99 




IQlO 


5000 


2.35 




109 


1000 


2.81 




109 


3000 


2.83 




IQlO 


500 


2.94 


Br?iup/Br7 


IQiO 


500 


0.66 




IQlO 


1000 


1.48 




iqs 


3000 


2.19 




iqs 


5000 


3.09 




IQlO 


3000 


3.34 




iqs 


1000 


3.45 


Br7/Panup 


IQlO 


1000 


0.44 




IQlO 


3000 


0.89 




IQiO 


500 


1.12 




109 


3000 


1.68 




109 


1000 


1.74 




IQlO 


5000 


2.04 
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determined from the scatter plots presented in Figure 2a. 
The model that best fits the data in Figure 4 (dashed 
Hue) has T = 1000 K and ne = 10^° cm'^. The pa- 
rameters and reduced x^ values for the best six models 
are presented in Table 3. We note that the tempera- 
ture sensitivity of these models at the best-fit density of 
10^*^ cm"'^ is most apparent at the mid-range ri„p-values 
where the most separation can be seen between the dif- 
ferent models. In contrast, the temperature is poorly 
constrained by the Pariup transitions with the lower and 
higher n,ip-values. 

Figure 5 contains the line ratio curve for the observed 
Brackett series features using Br7 as the reference tran- 
sition. Overplotted are four Case B model fits, all with a 
density of nc — 10^° cm^"^, which illustrate the temper- 
ature sensitivity of the models at the lowest n„p-values. 
The lowest temperature curves for T = 500 K and 1000 K 
are the two models that most closely fit the data. The re- 
duced x^-values for the six best fit models are presented 
in Table 3. Despite the rather large scatter bars associ- 
ated with these data points, the data significantly rule 
out the highest temperature Case B models. For both 
the Paschen and Brackett series line ratios, we find that 



Fig. 3. — Using the weighted mean line ratios for the Paschen 
series with Pa/3 as the reference transition (see Figure 2a), we plot 
the observed Paschen decrement curve as a solid line connecting 
the diamond-shaped data points. 'Scatter bars,' determined from 
the estimated variance, are plotted atop the data points. Of the 
possible 130 computed Case B models, we overplot the four with 
the smallest reduced x^-values. The temperatures and electron 
densities for these models are given in the figure legend. Reduced 
X^-values can be found in Table 3. 



In Figure 4, we plot a line ratio curve for the observed 
Paschen series features using Pa/? as the reference transi- 
tion and overplot the four Case B models with the small- 
est reduced x^-values. The 'scatter bars^' on each of the 
data points represent the estimated variance (see ^3.1^ 



We use the term 'scatter bars', instead of error bars, to more 



and T = 1000 K and 500 K for Paschen and Brackett, 
respectively. The consistency of the gas temperature and 
electron density ranges determined by both independent 
H 1 series line ratios and the statistical significance as- 
signed to these fits is remarkable and is strong evidence 
that the assumptions of the Case B approximation are 
valid for the emitting H i gas. 

Finally, it is instructive to compare the Case B model 
predictions to the observed Br7 to Pan^j, because of the 
large wavelength range covered by this ratio. In Fig- 
ure 6, we plot the observed line ratio curve for Br7 to 
Pan„p and overplot the best model fit with T = lO'^ K 

accurately reflect the nature of these these vertical bars. The range 
of values allowed by each 'scatter bar' is not due only to uncertainty 
in measured values, but also to the actual scatter of the line ratios 
(i.e., estimated variance) determined for all sources in our sample. 
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Fig. 4. — Using the weighted mean line ratios for the Brack- 
ett series with Br7 as the reference transition, we plot the ob- 
served Brackett decrement curve as a solid line connecting the 
diamond-shaped data points (see Figure 2b). Scatter bars, de- 
termined from the estimated variance are plotted atop the data 
points. Wc ovcrplot four Case B models of constant electron den- 
sity (ric =10^° cm~3). Note that the T = 500 K (dashed) and 
1000 K (dotted) models arc the two best-fit models (see reduced 
X^-values in Table 3). 
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Fig. 5. — The values of the line ratios of Br7 to the Panup of the 
Paschen series are plotted. The broad wavelength coverage between 
the Brackett and Paschen lines provides a stronger constraint on 
density as the predicted line ratio curves diverge from the observed 
curves significantly at high n„p. This sensitivity is demonstrated 
by plotting four models of constant temperature (T = 1000 K) 
and separated by one order of magnitude in density along with the 



data. The model with the smallest reduced x^ has Uc 
(see Table 3). 



10^ 



and no = 10 cm^'^. Additionally, four Case B line ratio 
curves, all with T = 10^ K (the temperature of the best 
model fit) and electron densities that bracket the best- 
fit density, no = 10^° cm""^, are also overplotted. At low 
n„p-values, the differences between the measured and the 
predicted line ratios are quite small for this temperature. 
However, for the line ratios with the largest Pariup states, 
the models with the same temperature and different den- 
sities are quite distinct from one another, clearly distin- 
guishing the model that most closely matches the data. 
As a result, the Br7 to Pan^p line ratios place a strong 
constraint on the n^ of the emitting gas. Once again, the 
temperature and density of the best-fit Case B model 
for Br7/Pan agree well with those determined for both 
the Paschen and Brackett series ratios. The large sepa- 
ration in wavelength between the Paschen features and 
Br7 permits this Case B analysis to be searched for any 
dependence on the reddening correction. The agreement 
between the model temperatures and electron densities 
determined from all three sets of line ratios suggests that 
the spectra have been properly corrected for reddening 
effects. 

3.2.3. x^ Contours 

In this section, the reduced x^-values for each of 
the 130 temperature and electron density dependent 
Case B models calculated for the Brackett series and the 
Br7/Pan line ratios are combined to determine which 
models are the best fits to the combined data^. In 
Figures 7a & 7b, the reduced x^ surface and contours 
are plotted, respectively, to illustrate the location of 
the models that are the best statistical fits to the data, 
amongst the grid of Case B models. In Figure 7a, the 
X^-values for every grid point are presented. The min- 
imum and surrounding surface locates the models that 
best fit the data. The steep rise of the surface indicates 
a sharp increase in the x^"Values, highlighting the nar- 
row range of temperatures and electron densities that fit 
the data well. Figure 7b is a contour plot of the confi- 
dence intervals associated with the models with the in- 
nermost contour representing the 60% confidence interval 
and the subsequent contours constrain the models out to 
the 99.9% confidence level. These confidence intervals 
also illustrate the tight constraints placed on both T and 
no. 

Figure 7 illustrates that the Case B analysis rules out 
an electron density much larger than 10^" cm~^, with 
10^^ cm~^ being completely ruled out. A density of 
< 10^ cm~^ is also ruled out at the 99% confidence level, 
for a temperature of 3000 K, and at an even higher con- 
fidence level for the rest of the temperatures considered. 
The models with densities of 10^° cm""^ have small re- 
duced x^ values, making this the most probable density 
for the emitting gas. For the gas temperature, a value 
greater than 5000 K is ruled out at the 99.9% level. It 
is clear that temperatures of < 2000 K are statistically 
favored, but temperatures lower than 500 K are also 
allowed and cannot be constrained without models for 
T < 500 K. A finer grid of models would be useful 
for better constraining the range of temperatures and 

^ Note that the x^ -values for the Paschen series model compar- 
isons were not included in this calculation to avoid introducing 
redundant data and artificially lowering the resulting x^-valucs. 
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Fig. 6. — a) The reduced x surface of the 13 density and tem- 
jerature dependent Case B models provided bv lStorev fc Hummer! 



Ill995al) are plotted. The valley of the surface locates the models 
with the smallest reduced x'^ and with its minimum approaching 
T = 1000 K and Uo = 10^" cm~^. b) The reduced x^ contours 
are plotted for the values of 1.049, 1.487, 1.666, 2.039, & 2.513 
corresponding to confidence intervals of 60, 90, 95, 99, & 99.9% on 
a T vs. Uc plot. The empty diamonds designate T, no grid points 
from the 18 Case B models which fall in the range of the plot. 
The contours illustrate the tight constraints placed on densities 
higher or much lower than 10 cm~^ and temperatures greater 
than 5000 K. 

electron densities predicted by the Case B comparison. 
Nonetheless, the analysis as presented here constrains 
both temperature and density to within a factor of a 
few. 

3.3. Comparison to Optically Thick and Thin LTE 

Cases 

In this section, we compare the observed line ratios to 
those predicted for gas at a variety of temperatures in 
both the optically thick (r > 1) or thin (t < 1) LTE 
regimes. 

3.3.1. Optically Thick Case 

In general, for optically thick gas in LTE, the intensi- 
ties of all lines approach the value for a blackbody of a 
single temperature. Assuming that the line widths are 
proportional to the wavelength of the line (as expected 
for Dopplcr broadening), the line flux ratios are deter- 
mined by 



fi _ Ai5x(Ai,r) 
F2 X2B\{X2,T) 



(1) 



where F is the line flux, B\{X, T) is the Planck function, 
A is the wavelength of the transition, and the subscripts 
1 and 2 distinguish between transitions in any ratio. 

Figure [7] compares the prediction of the optically thick 
LTE case to the observed Br7 to Pan^p line ratios. For 
the optically thick cases, a blackbody temperature range 
3000 < T < IC* K represents the range of expected 
temperatures in the accretion flow (e.g., Martin 199G; 
iMuzerolle et al.|[2001[ ). Figure [7] shows that the observed 
line ratios are not well-matched by the optically thick 
line ratio curves. In Table 4, the reduced x^-values cal- 
culated for several temperatures in the optically thick 
case, including those plotted in Figure [71 are presented. 
A temperature of w 3800 K is the best statistical fit to the 
data and has a reduced x^ of 5.4, placing it far outside 
the 99.9% confidence level. A similar result is found for 
the Paschen and Brackett series line ratios (not shown) 
where the minimum reduced x^ values are similarly ruled 
out. Therefore, this comparison demonstrates that the 
optically thick LTE case is not a good approximation of 
the H I emitting gas. 
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Fig. 7. — The Br7/Pa7i line ratio curves for gas in both the 
optically thick and thin LTE regimes are plotted. A comparison 
between the observed ratios and those expected for a gas in LTE 
demonstrates that the emitting gas is not in LTE. Shown are LTE 
cases for a variety of temperatures (as indicated in the inset), and 
there is no temperature that provides a good fit for either LTE 



case. The fact that the observed ratio at n^ 



■ 7 well-matches the 



LTE, optically thin value is direct evidence that the two observed 
lines in this ratio are optically thin. 





TABLE 4 




Optically Thicf 


AND 


Thin Reduced x^- 


VALUES 


Case 


T Reduced 




(K) 


2 

X 


Thick 


1000 


5e6 


Thick 


3000 


14 


Thick 


3800 


5.4 


Thick 


10000 


14 


Thin 


1000 


2400 


Thin 


10000 


32 


Thin 


100000 


18 


Thin 


CO 


17 



3.3.2. Optically Thin Case 

We next consider the optically thin LTE case, where 
the line flux ratios only depend on the relative popu- 
lation densities of the atoms in the upper n-levels, the 
associated Einstein A-values, and the wavelengths of 
the transitions. For a gas in LTE, the level popula- 
tions follow a Boltzmann distribution, which was used 
to calculate the optically thin LTE ratios for this part 
of the analysis. In Figure [71 the optically thin LTE 
curves for the Y^v^ /Variup line ratios with temperatures 
1000 < T < 10^ K are plotted for comparison to the 
curve representing the observed line ratios. 

As in the optically thick case, the observed data are 
not fit well by the line ratio values in the optically thin 
case. In order to quantify the quality of the fits, the 
reduced x^-values for the optically thin regime were cal- 
culated and are presented in Table 4. For comparison, 
the models with temperatures of 1000 K, lO** K, and 
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T^ cxo have reduced x^-values of ^ 10"^, 32, and 17, 
respectively. Once again, the LTE case falls outside of 
the 99.9% confidence level for any temperature. Simi- 
larly poor correlations for the independent Paschen and 
Brackett series line ratios also are found for the optically 
thin case. Therefore, this comparison demonstrates that 
the emitting gas is not in the optically thin LTE regime. 
By comparing the line ratio curves for the two high- 
est electron density Case B models (ug = 10^^ and 
10^^ cm~^) in Figure 6 to the optically thin LTE line 
ratio curves in Figure 8, it is clear that the Case B values 
approach those for an optically thin, high density gas in 
LTE. The convergence of these two regimes demonstrates 
how particle collisions dominate the statistics of n, Z-level 
populations in the Case B approximation at high den- 
sity. Thus, it is apparent that the emitting H i gas in the 
T Tauri systems in this sample is in a regime where colli- 
sions are important for the n„p-level populations, but not 
sufficient to give rise to a Boltzmann distribution. Only 
a sophisticated treatment of the level populations, such 
as that found in the current Case B theory, will properly 
predict the line ratios. 

3.3.3. Bry/PaS Ratio 

Upon inspection of Figures 6 and 8, we note that the 
Case B models and the optically thin LTE case at dif- 
ferent temperatures, predict nearly the same value for 
the Br7/Pa(5 line ratio. This ratio is unique amongst the 
others in this series because the electrons share the same 
n„p-level (n„j, — 7) . Thus the population density of the 
upper level of both Br7 and Pa6 is identical, and so this 
line ratio is independent of the manner in which the up- 
per state is (de)populated and is entirely independent of 
the temperature and density of the gas^. In the opti- 
cally thin limit, this line ratio is well-described by the 
following 



F2 



AiA 



1A2 



A2\l 



(2) 



where F represents the line flux, A is the Einstein A- 
value, and A is the wavelength of the transition. Only 
optical depth effects can cause the observed ratio to de- 
viate from that predicted by Equation [21 Therefore, the 
agreement between the observed value of the Brj/PaS 
line ratio and those predicted by Case B and the opti- 
cally thin regime is further evidence that the H i emission 
originates in an optically thin gas. 

At the same time, the value of the Br7/Pa(5 line ra- 
tio is accurately determined from a couple of well-known 
constants and, therefore, serves as a definitive test of our 
reddening correction and Case B analysis. An error in 
either the measurement of Av, the reddening law, or a 
systematic error in the determination of the line fluxes 
would lead to significantly different values measured for 
any of the line ratios. Therefore, the agreement between 
the predicted and measured values of the Br7/Pa^ is a 
verification of the data analysis techniques, the visual 
extinction values adopted for these sources, and the red- 
dening law employed. 

'' Note that the individual angular momentum states, I, asso- 
ciated with rtup- state are statistica lly populated for the densities 
considered here lIBrocklehurstllTQTlh 



4. DISCUSSION 

The Case B determination of gas temperature and elec- 
tron density of the emitting H i gas is completely inde- 
pendent of any assumptions about where the emission 
arises (e.g., in a wind, magnetically-channeled accretion 
flow, or disk). In this section, we discuss the T and Ue 
ranges constrained by the Case B analysis in the con- 
text of the magnetospheric accretion model, which has 
had the best success so far in explaining the origin of 
the H I emission features. In this picture, the emitting 
gas is confined primarily to the accretion funnel and is 
essentially freely falling onto the star. 

4.1. Accretion Flow Density 

For a fiducial accretion rate of M — lO^^A/© yr^^, 
magnetospheric accretion models predict a hydrogen 
density in the accretion flow to be ~ lO^'^ — 10^^ cm~'^ 
(|Martinlll996tlMuzerolle et al.lll998all200lD. For the ion- 
ization fractions calculated bv lMartiiini996f) . ~ 10~^ — 
10^"^, the range of electron densities for the accreting 
gas is expected to be 10^-10^^ cm"'^. As shown in sec- 
tion 13.2.31 the range of electron densities favored by the 
Case B analysis falls squarely within the range predicted 
by the magnetospheric accretion models. Therefore, the 
electron densities favored by the Case B analysis appear 
to be consistent with the expectations of the magneto- 
spheric accretion models. 

However, the tightly constrained values of the electron 
densities returned from the Case B analysis is represen- 
tative of 15 different stars possessing a range of spectral 
types and mass accretion rates. Given the large range 
in the mass accretion rates for the stars in this survey 
(see Table 1), which span more than two orders of mag- 
nitudes, it is surprising that the range of statistically 
significant electron densities in the Case B analysis is at 
most, a factor of a few p.2.3p . 

The fact that the large range of mass accretion rates 
is not reflected in the range of electron densities appears 
to require systematic source-to-source differences in the 
accretion flows, which serve to regulate Uc of the emitting 
H I gas. In a steady-state accretion flow, the hydrogen 
density is proportional to the mass accretion rate divided 
by the infall velocity and divided by the cross-sectional 
area of the flow. In the magnetospheric accretion models, 
the infall velocity approximately equals the escape speed 
from the stellar surface, and this is unlikely to span more 
than a factor of a few for the sources in this survey. Thus, 
there remain two potential processes for systematically 
regulating the electron density: 1) variations in the cross- 
sectional areas of the accreting flows and 2) fluctuations 
in the ionization fraction. 

Meas urements of the surface are a covered by accretion 
shocks (jCalvet fc Gullbrin3ll998( ) have the potential to 
be reliable probes of the cr oss-sectional area of the flows. 
ICalvet fcG uUbring ()1998f ) determined the areas of the 
accretion shock regions on the surfaces of the target stars 
in this survey to have a range of more than two orders of 
magnitude. Thus, the measured range in the area of the 
shocked regions, if representative of the variations in the 
cross-sectional areas of the flows, is sufficient to keep the 
electron density approximately constant for values of the 
mass accretion rates, which vary by a few orders of mag- 
nitude. Furthermore, the results of ICalvet fc Gullbrind 
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(|1998f) and lValentiet all (|1993l and C. Johns-KruU, pri- 
vate communication) are consistent with the area of the 
shocked regions being approximately proportional to M , 
a requirement for variations in the cross-sectional area of 
the flow to regulate the density of the accreting gas. Al- 
though the suggestion that variations in accretion rates 
are accompanied by nearly equal variations in the cross- 
sectional areas of the accreting flows is consistent with 
our finding of relatively constant Uc, this remains specu- 
lative. 

In addition to varying the cross-sectional area of the 
flow as a means of self-regulating the electron density, 
fluctuations in the ionization rate due to changes in the 
mass accretion rate and neutral gas density may also 
be responsible for maintaining a nearly constant electron 
density. Differences in M and the gas density may impact 
the recombination timescale and the ionization faction 
in such a way as to approximately regulate the electron 
density in the accreting flow. 

Of course, it is unclear whether either of the above 
processes are responsible for the small range of electron 
densities favored by the Case B analysis of the sample 
of TTS in this survey. Additional observational studies 
and further refinements of the line emission models will 
be necessary to address this issue. 

4.2. Accretion Flow Temperature and Luminosity 

The preferred values of the gas temperatures derived 
from the Case B analysis fall in the range of T < 2000 K 
and are substantially lower than the range of tempera- 
tures predicted by the magnetospheric accretio n mod- 
els o f 6000 5, T f 12000 K (see gU i MuzeroUe et al 
120011: iMartlij 119961 1 In the models of iMuzerolle et al 



()2001h . the gas temperature was a tunable parameter. 
The "optimal" temperatures of ~ 10"* K were found, such 
that the models reproduced the observed magnitudes of 
the H I emission line fluxes. Martin (1996) presented 
a self-consistent calculation of the thermal structure of 
the accretion flow, showing that adiabatic compression 
was the strongest heat source considered for the accret- 
ing gas. That work found that the heating and cool- 
ing rates reached equilibrium for a gas temperature of 
T ^ 6000 K, and the total luminosity of the accre- 
tion flow is on the order of 10^^-10'^" ergs s^^. This 
value is comparable to the line luminosities of some of 
the strongest H i emission features (e.g., Pa/3, Vaq, Br7; 
Figure [2) observed for sources in this survey (assuming 
D ~ 140 pc) and cannot account for the total luminosity 
of the accretion flow. While a higher gas temperature 
of ^ 10"* K provides a better estimate of the total lu- 
minosity of the flow, it requires more heating than that 
provided by the heat sources included in Martin (1996'). 
Therefore, the primary mechanism responsible for pow- 
ering the luminosity of the accretion flows remains un- 
clear. In addition, the low gas temperatures returned by 
the Case B analysis raises questions about the cooling 
rate of the flow. 

4.2.1. Cooling Time of Gas in the Accretion Flow 

In order to gain more insight into the thermal struc- 
ture of the accretion flow, in this section, we estimate the 
radiative cooling time of the gas, TcooI- IMuzerolle et al.l 
(|2001f ) found that the total radiated luminosity of the 



accretion flow, Laf, is comprised primarily of hydrogen 
continuum emission (> 75%), plus a significant contri- 
bution from emission line fluxes. In a comparison of iaf 
to the accretion luminosity, Lace, assuming that 



Faf — jFg^cc — J 



GAFM 



(3) 



where M* and i?* are the stellar mass and radius, and / 
is the fraction of the accretion lumino s ity ra diated away 
by the accreting gas, IMuzerolle et al.l (|2001h found / to 
be ~ 5%-30%, for the range of accretion rates exhibited 
by our sample stars. 

The total mass of material in the accretion flow at any 
instant in time can be approximated by Mrg, where rg 
is the time for material to free-fall from the inner disk 
truncation radius, Rt, to the stellar surface. Tff is well 

approximated by R^ (GM*)^^/^, which is several hours 
for a typical cTTS. The total thermal energy of the gas 
in the accretion flow is 



3 MTff 

Uth '^ -zi^T , 

2 /ifrifj 



(4) 



where k is the Boltzmann constant, ^ is the mean molec- 
ular weight, and ttih is the proton mass. 

The cooling time, TcooI, is given by, f/th/-^af- It is in- 
structive to compare the duration of the cooling time 
with that of the free-fall time by using the ratio 



T-cooi nnn-^/^0.lW^_W200kn^\ 

where the input parameters assume fi — 1. For the 
fiducial values in this equation, the cooling time is much 
shorter than the dynamical (free-fall) time. 

This appears to be true for most TTS, since Equation 
© is nearly independent of M and Rt (the onl y depe n- 
dence is via the parameter /; IMuzerolle et al.l (j2001h ). 
This result is the consequence of the gas having a lumi- 
nosity equal to a substantial fraction (/) of the gravi- 
tational potential energy released (i.e., the accretion lu- 
minosity), while having a temperature much cooler than 
the virial temperature. Therefore, if the accretion flow 
is indeed radiating with a luminosity such that / ^ 0.1, 
the free-falling gas has time to reach an equilibrium tem- 
perature, set by a balance between heating and radiative 
cooling. This will be true, even for the higher tempera- 
tures predicted by magnetospheric accretion models. 

In order to simultaneously explain the large observed 
luminosities and low temperatures favored by the Case B 
analysis, the accretion flow radiation models need to 
be improved. In order to power the emissio n, addi- 
tional heating is required than was considered by lMartinl 
(1996). Furthermore, the radiative cooling of the gas 
must be more efficient than was found by any of the pre- 
vious models, to explain how the observed line fluxes can 
be produced by a gas with a lower equilibrium tempera- 
ture. 

4.2.2. How Important is Ionizing Radiation? 
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The models of iMuzerolle et all (|2001h and iMartinI 
()1996D included an ionizing source of photons with a lu- 
minosity of '^ 10^^ erg s"'^ (shortward of 912 A). Since 
this is not capable of powering the observed line lumi- 
nosities, it is not surprising that these works favored 
temperatures on the order of ~ 10 K, where thermal 
(collisional) excita tion/ionization effects become impor- 
tant for hydrogen (jMuzeroUe et al.lll998af ). However, the 
temperature favored by the Case B analysis is well-below 
this temperature, suggesting that a non-thermal ioniza- 
tion mechanism is at work. 

High energy photons generated in both the active coro- 
nae and the accretion shocks at the surfaces of active 
TTS are potential sources of incident radiation capable 
of ionizing the accreting gas. The observed X-ray lu- 
minosities of the stars in this sample are in the range 
10^^ — 10'^^ erg s~^ (see Table [T]). The luminosities in 
the extreme UV (shortward of 912 A) may be similar in 
strength or potent ially even higher (Herczeg et al. 2002; 
iBergin et al.ll200^ . Furthermore, depending on the ge- 
ometry of the accretion columns and the mass accretion 
rates, the hydrogen column densities of the accretion 
flows are estimated to range between ~ 10^" — 10^^ cm~^. 
At these column densities, the accreting gas will be capa- 
ble of absorbing both UV and X-ray photons emitted in 
the nearby coronae or from accretion shocks on the stel- 
lar surface. The absorption of this radiation, which may 
be a significant fraction of the high energy photons pro- 
duced (Gregory et al. 2007) has the potential to power 
much of the H i emission. 

A higher incident ionizing flux than that included in 
previous models should predict a higher ionization frac- 
tion for a given temperature. In turn, this will lead to an 
increase in the cooling rates from e.g., free- free emission 
and radiative recombination. Thus, it qualitatively ap- 
pears that an irradiated accretion flow model including a 
realistic ionizing flux may favor a temperature well-below 
~ 10^ K and be more consistent with the temperature 
indicated by the Case B analysis of this work. 

5. SUMMARY AND CONCLUSIONS 

In this multi-epoch spectroscopic survey, the line ratios 
from 17 Paschen and Brackett H i emission features were 
measured from a total of 73 spectra of 15 actively accret- 
ing cTTS. The values of the line fluxes measured for the 
strongest features (e.g., Pa/3) over the entire sample span 
almost two orders of magnitude. In spite of this large 
variation in line fluxes, the measured values of all line 
ratios were approximately constant, within the observa- 
tional uncertainties, from source to source and epoch to 
epoch. These results are consistent with earlier work that 
also found little variation in the value of the Pa/3/Br7 ra- 
tio (jMuzerolle et al.ll200lHNatta et al.ll2006f) . Due to the 
little variation observed in the values of the line ratios, a 
single value, given by the uncertainty-weighted average 
of all the data for each ratio, was adopted. 

5.1. Level Populations and Optical Depth 

The measured values for the H i line ratios and those 
expected for both the optically thick and thin LTE cases 
for any temperature were not in agreement. This result 
is evidence that the n, Z-level populations of the H i gas 
are not in LTE. On the other hand, the observed ra- 
tios were statistically well-fit by the predictions of the 



recombination line theory Case B approximation, tightly 
constraining the values of T and Uq. This is evidence 
that the characteristic assumptions of Case B theory hold 
for the emitting gas. In addition to the success of the 
Case B models, the measured value of the temperature- 
and density-independent ratio of Br7/Pa^ is consistent 
with the optically thin value (i j3.3.3p . 

5.2. Density 

The observed line ratios are statistically well-fit 
by Case B models with an electron density of 
nc = 10^° cm"'^ fi )3.2.3[) . Furthermore, electron densi- 
ties of > 10^^ cm~'^ and < 10^ cm~^ are statistically 
ruled out at the 99% level or higher (depending on the 
temperature) . Thus it appears that the range of electron 
densities allowed for all 15 stars in this survey and for all 
epochs is constrained to lie within only a factor of a few 
of 10^° cm-3. 

The electron density favored by the Case B analy- 
sis is consistent w ith that pred i cted by magnetospheri c 
accretion models ()Martinlll99^ IMuzerolle et al.lll998a[ ). 
for typical cTTS parameters and assuming an ionization 
fraction of ~ 10~^ — 10~^ (SH)). However, it is a sur- 
prise that the allowed range in electron densities is so 
narrow, given that the range of accretion rates exhibited 
by the stars in our sample spans more than a factor of 
100. If the H I emission arises primarily from material 
in a magnetospheric accretion flow, this seems to require 
systematic differences from source to source that regulate 
the electron density (see Q. 

5.3. Temperature 

The Case B statistical analysis rules out electron tem- 
peratures of T > 5000 K for the emitting H i gas, at the 
99.9% level fJ lS^S]) . Temperatures of < 2000 K are fa- 
vored, and temperatures as low as a few hundred Kelvin 
cannot be ruled out. These temperatures are well below 
the predictions of the magnetospheric accretion models 
of 60 00 < T < 12000 K (MuzeroUe et al. 2001 : Marti^ 
Il996[ ). However, these models currently do not include 
all important sources of heating, cooling, and ionization 
( §4.2p . Therefore, the lower gas temperatures we find 
may still be consistent with H i arising from accreting 
gas. 

Independent of the location of the H I gas, the Case B 
temperature range appears to be too low to produce the 
observed emission line luminosities, unless a non-thermal 
process is ionizing the gas f ii4.2.2p . We propose that 
a significant flux of ionizing photons, arising from the 
stellar coronae and/or accretion shocks near the stellar 
surfaces, may simultaneously explain the high observed 
line luminosities and low temperatures of the emitting 
H I gas. 
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